Spatial entanglement of bosons in thermal equilibrium 
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We investigate spatial entanglement in a system of spinless non-interacting bosons in a one di- 
mensional harmonic trap described by the grand canonical ensemble. We find the lower bound for 
the amount of spatial entanglement contained in the average state of two bosons in the trap and 
show that it exists for arbitrary high temperatures. 
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\Q . I. INTRODUCTION 

o : 

Since its discovery, the phenomenon of entanglement lacks full understanding. For instance, when one deals with 
mixed states, we do not have necessary and sufficient criteria to determine if a given state is entangled or separable. 
One has, of course, numerous criteria for entanglement, for example the partial transposition test or the reduction 
criterion 0, but they provide only necessary conditions. Mixed states, on the other hand, are commonly observed, as 
quantum systems interact with their environment. In particular, a natural state of the matter is a thermal mixture. 
y—{ ' Therefore, the question whether a (thermal) density matrix is entangled or not, and how entanglement is related to 
£NJ , thermodynamical principles and quantities can arise in many physical situations. The notion of entanglement becomes 
even less clear when the effects of Bose-Einstein (BE) or Fermi-Dirac (FD) statistics cannot be neglected. This is 
the case when the wave functions of particles forminga physical system significantly overlap (as in thermal mixture). 
This problem has recently been investigated in Ref. 0,|j,l3)0)0)lEI3|- 

There are many interesting physical situations where BE and FD statistics play an important role. Standard 
examples are fermionic and bosonic gases, quantum phase transitions |Toj . super-conductivity super-fluidity 
j— ( ' ^2 an d Bose-Einstein Condensation (BEC) |13|. a phenomenon in which a large group of indistinguishable bosons 
is described by the same single-boson ground-state wave function. It is interesting to see whether entanglement is 
related to these phenomena and to what extent. Concerning BEC, this problem was addressed for zero temperature 
in the Ref. 0] and for non-zero temperature in the Ref. ^3 . 

In this paper we discuss the so-called spatial entanglement in the simplest model of non-interacting spinless 
f> |] bosons trapped in a one-dimensional harmonic oscillator potential described by the grand canonical ensemble. In 
i 1 > this example BE statistics cannot be neglected and the system exhibits phase transition in the form of Bose-Einstein 
£h ' condensation |Tt3 |. 

=5 

II. SPATIAL ENTANGLEMENT 

> ' 
• i— I , 

We start our discussion by defining what we understand by spatial entanglement in a system of indistinguishable 
non-interacting bosons. We limit our discussion to a one dimensional case but the extension to two and three 
?3 ! dimensions is straightforward. 

Let us first consider the simples possible case of two bosons moving in some trapping potential. The eigenkets of 
the Hamiltonian corresponding to eigenvalues Ek are given by \<pk) (k = 0, 1, 2, . . .) and the symmetric basis for two 
bosons reads \ip)k,i — 1 ^ fci (\4>k) \4>i) + \<Pi) \4>k)) + o~ki\(t>k) \4>k)- Are these basis states entangled? This question does not 
make sense unless we precisely define what the degrees of freedom are in relation to which we discuss entanglement. 
Additionally, bosons are indistinguishable, which further complicates the answer. 

To see the problem more clearly, let us consider a state describing two bosons that are spatially separated, i.e., their 
respective wave functions do not overlap. A simple example is given by two indistinguishable bosons each confined 
in its own infinite square well potential, say A and B. The cigenstates of the wells A and B are \<jr£) and \<j>f ) 
respectively. It is possible to prepare the particles in the wells in the state -L(|a^)|0f ) + l^f)]^)) with k ^ I. This 
state is, without doubt, entangled since the the partially transposed density matrix with respect to the subsystem 
B has negative eigenvalues. Here there is no ambiguity which degrees of freedom (subsystems) are entanglement or 
in other words what the subsystems A and B are. We can unambiguously identify such subsystems because the two 
bosons are spatially separated and become physically distinguishable in the sense that local measurements performed 
on them will not reveal the bosonic properties of the particles. 

Now, let us put these two bosons into the same trapping potential such that their energies are Ek and Ei with 
k 7^ I. The state describing this situation is ^(|</>fc)|</>z) + \(j>i)\<j>k)) an d it looks deceptively similar to the previous 
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one. However, now there is nothing that makes these two bosons different or, as opposed to the previous example, 
there are no well-defined subsystems A and B. We now investigate what happens to entanglement in this case. 

We could naively define subsystems A and B by analogy with the previous example, i.e., the first ket would 
correspond to A and the second one to B. Then we could perform partial transposition with respect to A and check 
for the negativity of the density matrix obtained this way. We would get exactly the same result as in the previous 
example. However, this operation is physically meaningless because the result we obtain is a negative density matrix 
that is not symmetric with respect to a relabeling of the particle indices. This is a serious problem because partial 
transposition corresponds to the reversal of time on one of the subsystems, which cannot change the fact that we deal 
with bosons. This in turn means that even the unphysical negative operator after partial transposition must be still 
symmetric. We did not have this problem with the previous example where the bosonic character of the particles was 
irrelevant due to the physically meaningful labels distinguishing them. The labels enumerating bosons in this example 
do not have any physical meaning (there is no boson one and two) and one has to find other means to distinguish 
them. We have mentioned this example in relation to the Ref. ^3 where the criterion of entanglement for two bosons 
in a harmonic trap was based on the partial transposition with respect to the indices enumerating bosons. In view of 
the above arguments the results in the Ref. 0] are not correct. 

In the following part of this section we define physically meaningful subsystems (first for two then for many bosons) 
A and B and discuss entanglement between them. In order to distinguish this type of entanglement from other types 
of entanglement we henceforth refer to it as spatial entanglement. Our definition is not unique and one could conceive 
other ways of defining entanglement |l8j . The reason why we choose to investigate the spatial entanglement defined in 
this way is that it may be related to the phenomenon of off-diagonal long range correlations, which in turn is related 
to phase transitions jl9j . 

First we define two non-overlaping adjacent intervals dividing the whole space into two parts, say A and B. With 
respect to this division the eigenket \<pk) in the position representation can be now written as (x\<f>k) = ^fe(^) = 
^/p~A{k)4 l k ( x ) + VPsity^k ( x ) i where ^Jpa(b) W) ^ s * ne probability of finding the particle in the region A and fit^ ( x ) 
is the wave function defined as 1 4>k(x)xA(B)(x). Here XA(B)( X ) is a characteristic function of the region A(B). 

i.e., it takes the value 1 if a; belongs to A(B) and otherwise. Similarly, the symmetrized state ipkl{ x ii x "i) °f two 
bosons (in position representation) can be written as 

ipki{xi,x 2 ) = \/pAA(kl)ipki A (xi, x 2 ) + VPBB(kl)ipki B (xi,x 2 ) + ^pAB(kl)ip k A i B (x 1 ,x 2 ), (1) 

where 

^kl A ( x ^ x ^) = i fl 1 ^ h{xx,x 2 )xa{xi)xa{x2) 

y/PAA{kl) 



tf>kl B ( x l' X 2) = i M 1 . 4>kl{xi,X 2 )XB{xi)xB{X2) 
\/PBB(kl) 

^ki B ( x i, x 2) = —^===iPm(xi,x 2 )(xa(xi)xb(x 2 )+xa(x 2 )xb(xi)). (2) 

VPAB{kl) 

Obviously, p A A(kl) denotes the probability of finding two bosons in A, pAB(kl) denotes the probability of finding 
one of the bosons in A and the other one in B etc. These probabilities can easily be found by computing the 
appropriate integrals. All the above wave functions are symmetric and they describe sub-ensembles corresponding 
to three physically distinguishable situations; two bosons in the region A, two bosons in the region B and bosons in 
separate regions. 

Let us now discuss the spatial entanglement between these two regions. We can compute the amount of entanglement 
contained in the state \ipki) by computing the entropy of the subsystem A, which results from tracing out the subsystem 
B (this is a good measure because we deal with pure bipartite states) . This is a similar approach to that in Ref. . 

Before computing the von Neumann entropy of the subsystem A let us write the state Q without using the position 
representation. The notation we adopt here is similar to that of the occupation number representation in the so-called 
second quantization. We consider the two cases k ^ I and k — I separately. For k ^ / we have 



|^w) = y/PAA(kl)\l k , h) A \0) B + y/pB B (kl)\0) A \l k , h) B + 

WpA{k)p B {l)\lk) A \h)B + y/PA(l)PB(k)\h) A \lk)B), (3) 

where |0),4(b) denotes the vacuum in the region A{B). This is not a physical vacuum like in field theory, it is simply 
a ket conveying the information that there are no particles in A(B). The state means that there are 

two bosons in the region A(B) one of them in the state l^t^) an( i the second one in the state {(pf^). The states 
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1/c)a|1z) b and |I;)a|1/c)b describe the situation when one of the bosons is in the region A and the second one in the 

A(B 

fe(i) 



region B in the appropriate states |</v/,f ^). More precisely 



|U,i«Ub)= 7 1 ,, n (l^ (B) )l^ (B) ) + l^ (B) )l^ (B) )) 

V-PAA(BB)(K<) 

|l^|l fe ) s = i=(|^)|^) + |^)|^)). (4) 

It is important to stress that the normalized states \l k , h)A(B) and |1zc)a(b) do not form an orthogonal basis in the 
two- and single-particle subspace respectively in the Hilbcrt space describing the sector A(B). We also note that these 
states are symmetric with respect to exchange of fictitious indices enumerating bosons. 
It is interesting to write explicitly the probabilities appearing in the above formulas 

p AA (kl)=p A (k) PA (l)(l + \ A (l k \h) A \ 2 ) 
p BB (kl) = p B {k)p B {l){\ + |B(lfe|l ( )B| 2 ) 

PAB(kl) = 1 - p AA {kl) - PBB(kl). (5) 

We observe a bunching effect since the probability of finding two bosons either on the left or on the right is enhanced 
in comparison to distinguishable particles for which the respective probabilities would be p AA (kl) = p A {k)p A {l) and 
p BB (kl) = ps(fc)ps(0- This enhancement means that the probability of anti-bunching, i.e., finding each boson 
in a separate region is reduced when compared to distinguishable particles. The effect disappears if the overlaps 
U(B)(lfc|l;)yt(B)| are zero. If k = I the formulas simplify considerably and we have 

\il>kk) =PA(k)\2 k ) A \0) B + PB (k)\0) A \2 k ) B + y/2 PA (k)p B (k)\l k )A\lk)B, (6) 

where 

|2 fe )^ ( B) = |^ (B) )l^ (S) ) 
1 
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l k )A\l k )B = ^=(\^)\^) + \^M)). (7) 



Again, the state \ipkk) as well as all its components are symmetric with respect to fictitious indices enumerating 
bosons. 

A remark concerning the notation used is in order here. The states |lfc)^i|li)s, ^fc^lO)^, 10)^12^.)^ or |lfc)A|lfc)s 
are separable with respect to the subsystems A and B in spite of their rather complicated form when written using 
original kets. This is because they are invariant under partial transposition with respect to the subsystem A or B. 
Here partial transposition, say with respect to the subsystem B, simply amounts to exchanging indices between kets 
and bras having label B. For instance, (\(f> k )\4>m)( ( t ) n\(4't\) TB = \4>k )l0n ){ ( t>m\{ ( t > f\- This ensures that after partial 
transposition the state remains symmetric with respect to exchange of indices enumerating bosons. 

We clearly see the tensor product structure of two subsystems A and B. After carrying out the trace over the 
subsystem B we get 

p A (kl) =PAA(kl)\l k ,h) AA (lk,h\ +p BB (kl)\0) A A(0\ +p A (k)p B (l)\l k ) AA (l k \+p A (l)p B (k)\h) AA (h\ + 
VPA(k)p B (l)p A (l)p B (k)(\l k ) AA (h\ + \h) AA (l k \) (8) 

for k ^ I and 

p A (kk) =p 2 A {k)\2 k ) AA (2 k \ +p 2 B (k)\0) AA (0\+2 PA (k)p B (k)\l k ) AA (l k \ (9) 

for k = I. 

Let us now discuss the consequences of the adopted definition for the amount of spatial entanglement. First of 
all, the states \ipkk) are entangled with the amount of entanglement given by the Shannon entropy of the binary 
probability distribution p A (k) 2 ,p B (k) 2 , 2p A (k)p B (k). Please note that according to the incorrect reasoning presented 
before when one identifies subsystems with bosons themselves, there is no entanglement contained in this state. 
Similarly, the states with different k and I are not maximally entangled unless the effect of bunching disappears. 

The amount of spatial entanglement between the regions A and B depends on the way the whole space available 
to the bosons is divided. For instance, if the bosons are trapped in a one dimensional harmonic trap a division of 
the whole space into two symmetric parts (with the demarcation line going through zero) yields a different amount 
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of entanglement then an asymetric division. In the example of a one dimensional harmonic oscillator the symmetric 
split yields the largest amount of entanglement. 

It is straightforward to extend the above discussion to many-boson systems. We are not going to discuss this case 
here as in this paper we always deal with two bosons in the form of a reduced density matrix of the bigger many-boson 
system. Finally, we would like to point out that one can also define spatial entanglement in the sense specified for 
non-interacting fcrmions. 

In the next part of the paper we study the properties of spatial entanglement between A and B as defined above 
in a system of non-interacting bosons in a grand canonical ensemble in a harmonic one dimensional trap. 



III. SPATIAL ENTANGLEMENT IN THE GRAND CANONICAL ENSEMBLE 



We consider a system of non-interacting spinless bosons trapped in a one dimensional harmonic oscillator potential 
in thermal contact with a large system (reservoir). If we allow for the exchange of bosons between the trap and the 
reservoir the density matrix describing bosons inside the trap is called the grand canonical ensemble and reads 



1 x x 

i ,= 7- E E ■•■ ex p 

no— ni— 



Kk=0 




\n ,ni,...)(n ,m,...\, (10) 



where T is the temperature of the reservoir, k is Boltzmann constant, E k the energies of harmonic oscillator, p 
chemical potential and the ket |no,ni, . . .) describes a symmetrized state of no bosons occupying the ground state, 
rii bosons occupying the first excited state etc. In all the calculations we put n = 1 and E k = k (fc = 0, 1, . . . ). 

Due to the exchange of bosons between the trap and the reservoir the number of bosons in the trap fluctuates. To 
establish a link between experiments where the number of bosons in the trap is constant (such a situation is described 
by the canonical ensemble) one has to impose a condition that the mean number of bosons in the trap (N) is constant. 
This can be done by adjusting chemical potential with temperature. Unfortunately, to find a required dependence of 
fi on temperature one has to resort to numerical methods. 

One can now apply the method presented in the previous section to discuss the spatial entanglement in the system. 
This is a very difficult task because of the large number of bosons in the trap (this number fluctuates around the mean 
value (N)). However, one can ask an interesting question of how much of spatial entanglement is contained in an 
average pair of bosons in the trap. The amount of this spatial entanglement may be related to the phenomenon of the 
so-called off-diagonal long-range correlations, which indicate the onset of Bose-Einstein condensation (BEC) (or, in 
general, phase transitions) |19|. In the system considered here (finite and one dimensional system of non- interacting 
bosons) one has to be careful how to define BEC because, strictly speaking, it does not exist. However, as discussed 
in Ref. |l6j| one can observe a phenomenon of the macroscopic occupation of the ground state, which may be taken 
as the definition of BEC. 

Therefore, we focus our attention on the reduced density matrix p(2) of an average pair of bosons in the trap 

fc k>l 

where M — \ (3 X)fc( n fc) 2 + (^) 2 ) > (N) — Yl k { n k) ^ s the mean number of particles in the trap and (rife) 
: exp[^(E k — p,)] — 1) _1 is the mean number of bosons occupying the energy eigenket \(j>k) of harmonic oscillator 
potential. 

To find the amount of the spatial entanglement contained in p{2) we could compute the so-called negativity in- 
troduced in Ref. |2fJ. The negativity is the sum of moduli of the negative eigenvalues of the partially transposed 
with respect to B (or A since the system is symmetric) density operator p(2). Unfortunately, it is very difficult to 
diagonalize p(2) Ts since it is an infinite matrix. One would have to cut off the summations in (| 1 1 1> and diagonalize it 
numerically but a reasonable cut-off leads to a huge matrix which dimension increases too fast to be tractable numer- 
ically. However, we can find an analytical formula for the lower bound for the total spatial entanglement contained in 
p(2). This lower bound is related to the spatial entanglement related to a "coherent bunching" of bosons in either A 
or B, i.e., to the states of the type 1i)a\0)b + |0)^| 1^ , 1i)b- To see this let us carry out the partial transposition 
of the matrix p(2) with respect to B 



P(2) TB = E^ (|2fcU|2fc) S A<0| B <0| + |0) j4 |0) BA <2 fc | B (2 fc |) + 
fc 

^-J2PAA(k,l)(n k )(m)[(\l k , h) A \lk, 1 ; )ba(0| b (0| + \0) A \0)B A (lk, h\ B (lk, h\) + A(2), (12) 



M 

k>l 
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where A(2) is orthogonal to the first two operators and is not, in general, a positive operator. Therefore, the negatives 
eigenvalues of the first two terms give us the lower bound for the total entanglement in p(2). Taking into account the 
possible negative eigenvalues of A(2) can only increase the amount of total spatial entanglement. 
To find the negative eigenvalues of the first two operators in l|12f> let us define a ket 



IX) = \(n k ) 2 \2 k ) A \2 k ) B + X)<»k) WIU, h)A\U, h) B - (13) 

With this notation we have 



2 

fe k>i 



p(2) Tb =^(|x)<0| + |0)(x|) + A(2), (14) 



where the ket |0) = |0)a|0) _b- The operator |x)(0| + |0)(x| has only one negative eigenvalue A given by 

(n k ) 2 (n l ) 2 A (l k \l l ) i A + 



\ k,l 

^ X] ( n k)( n l)( n l')( n k')A(^k\^k'}AA(h\h')APAA(kl)pAA(k'l') + 
k>l k'>l' 

X)< n *) 2 < n *'X n '')A< 1 fcl 1 *'>A>l(lfc|li'>APAA(* / 0) , (15) 



k'>V k 



where we have used the symmetry of the division into A and B. Therefore, the lower bound for the negativity of the 
state p(2) equals to A. Its behaviour with temperature is depicted in FIG. 1. 




FIG. 1: Plot of the lower bound for the negativity of the state p(2) (given by A) versus temperature for the fixed mean value 
of the particles (N) . 



IV. CONCLUSIONS 



From the formula for A it can be easily shown that the spatial entanglement of an average pair of two bosons in 
the trap exists for arbitrary temperatures and it does not vanish with the increasing (N) . This strongly suggests (see 
discussion in the Ref. [2l]]) that the multipartite spatial entanglement of the bosons in the trap, i.e., the total spatial 
entanglement contained in the state p, should exist for arbitrary high temperatures. 

It would be interesting to compute the total spatial entanglement contained in the state p(2) and investigate its 
behaviour around the critical temperature for Bose-Einstein condensation (BEC) understood as the macroscopic 
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occupation of the ground state (see a detailed discussion in the Ref. 16]). We conjecture that there should be a 
relatively sharp decrease of the amount of spatial entanglement contained in p(2) around this temperature. This, in 
turn, would suggest that there should be an analogous sharp change of the total spatial entanglement contained in p. 

As a final remark we would like to point out that entanglement related to the "coherent bunching" of bosons in 
A and B, i.e., the spatial entanglement contained in the states 1i)a\0)b + |0)yi|lfc, h)s is similar to the so-called 
non-locality of a single photon [22] ■ 
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